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ABSTRACT 

A subgrid-scale model recently derived by Yoshizawa* for use in the 
large-eddy simulation of compressible turbulent flows is examined from a fund- 
amental theoretical and computational standpoint. It is demonstrated that 
this model, which is only applicable to compressible turbulent flows in the 
limit of small density fluctuations, correlates somewhat poorly with the re- 
sults of direct numerical simulations of compressible isotropic turbulence at 
low Mach numbers. An alternative model, based on Favre-filtered fields, is 
suggested which appears to reduce these limitations. 
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In a recent article, Yoshizawa* developed sugrid-scale models for 
possible future use in the large-eddy simulation of compressible turbulent 
flows. These models were obtained by making use of a multiscale version of 
Kraichnan's DIA formalism. 3 While the study conducted by Yoshizawa is quite 
interesting, it will be pointed out in this work that the specific subgrid- 
scale stress model derived therein is limited to slightly compressible flows 
(i.e., to flows with small density fluctuations) and also correlates somewhat 
poorly with direct numerical simulations of isotropic compressible turbulence 
at low Mach numbers. This poor correlation arises from Yoshizawa's formula- 
tion of a Smagorinsky-type model which neglects entirely the effects of momen- 
tum exchanges between the small and large scales that are accounted for in the 
subgrid-scale Leonard and Cross stress terms. These inaccuracies can be re- 
duced by a new subgrid-scale stress model that was recently derived by the 
authors.^ The purpose of this work is to provide a comparison of these two 
alternative approaches. Subgrid-scale models are needed for the large-eddy 
simulation of compressible turbulence with important applications to high 
speed aerodynamic flows. 

In the large-eddy simulation of compressible turbulent flows, any field 
quantity F in the flow domain D can be decomposed as follows 

F « F + F" (1) 

where 

F - / G(x - z,A)F(z)d 3 z (2) 

' /V /V /w 

D 

is the filtered part and T" is the subgrid scale part which accounts for 
the scales that are not resolved by the computational grid A* Yoshizawa* 



developed subgrid scale models (based on the filtering procedure (1)) by mak- 
ing use of Kraichnan's DIA formalism generalized to the compressible regime 
with a multiscale expansion. It was assumed therein that the length and time 
scales of the fluctuating fields were small compared with those of the mean 
fields. Such an assumption restricts the applicability of the models to flows 
with small density fluctuations (i.e., it is equivalent to making a small com- 
pressibility expansion). A Smagorinsky model for the subgrid scale stress was 
obtained to the first order when a state of isotropy was assumed for the zero- 
order solution. Comparable subgrid scale models of the gradient transport 
type were derived, to the first order, for the remaining subgrid scale corre- 
lations (e.g., the subgrid scale heat flux). 

As an alternative to this approach, in line with the more traditional 
analyses of compressible turbulent flows, a Favre filter can be defined as 
follows : 

?-? (3) 

p 

where p is the mass density. This gives rise to the alternative decompo- 
sition 

F - F + F' (4) 


where F' is the subgrid scale part of F based on Favre filtering. The 

Favre filtered equations of motion for an ideal gas are given by^ 
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where v k is the velocity vector, p = pRT is the thermodynamic pressure, 
T is the temperature, and R is the ideal gas constant. In (6) - (7), 


2 3v k 3v Jt 

V " - 3 "5-Z 6 k£ . + + F^> 


( 8 ) 


T kt " - t( Y« - \\ + 'A + V* + v k v p 
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Q k - c p p( V - V + v k * + V' + v k T '> 


(11) 


are the molecular stress tensor, the subgrid scale stress tensor, the viscous 
dissipation, and the subgrid scale heat flux, respectively, whereas y is 
the dynamic viscosity and < is the thermal conductivity. 

In order to achieve closure of the equations of motion (5) - (7), subgrid 
scale models for and Q k must be provided. We recently proposed the 
following models^ 
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(13) 


(14) 


1Z ~ - s u s u- <l5 > 

S 

Pr-j. is the turbulent Prandtl number (which assumes a value of approximately 
1/2), and C R and Cj are dimensionless constants which were found to assume 
the values of 0,012 and 0,0066, respectively, by correlating with the results 
of direct numerical simulations of compressible isotropic turbulence. This 
model may be viewed as a compressible analogue of the linear combination model 
of Bardina, et al.^. The subgrid scale stress model (12) was found to corre- 
late well with the results of direct numerical simulations of compressible 
isotropic turbulence for average Mach numbers 0 Mq _< 0.6 (for average 
Mach numbers there can be localized supersonic regions in this 

flow^). Correlation coefficients in the range of 80-95% were obtained at the 
tensor level independent of the average Mach number (see Table 1). These 
direct simulations of isotropic turbulence were performed using Fourier collo- 

O 

cation methods on a 96 J grid at a Reynolds number of Re^ = 40 based on 
the Taylor microscale. The initial energy spectrum was chosen to match that 
measured by Comte-Bellot and Corrsin^, and the results (including those to be 
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shown later in Fig. 1) were then non-dimensionalized with respect to a 

reference length scale and velocity scale of 10/ir cm and 1 cm/sec, 

respectively (the physical properties of air at room temperature were used for 

the reference density and viscosity) . The rms density fluctuations 

p = <p' 2 > 1/2 /<p> for M n - 0.1, M n - 0.4, and M n - 0.6 were 0.0034, 
rms u u u 

0.0358, and 0.0739 respectively, where <•> denotes the spatial average. 

In contrast to these compressible subgrid-scale models based on Favre- 
filtering, Yoshizawa* represents the subgrid scale Reynolds stress v''v'' 
by the Smagorinsky model (derived formally to the first-order using a two- 
scale DIA method; see equation (81) of Yoshizawa*) 
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and ^ uu 2 is a dimensionless constant that was found to assume a value of 
0.16. However, by taking the trace of (16) it follows that 

V*v - 0 (19) 

fS* /w 


which is strictly true only for incompressible turbulence. Hence, as alluded 
to earlier, the subgrid scale model (16) of Yoshizawa 1 is only appropriate for 
weakly compressible turbulent flows (this limitation arises since 


- 6 - 


a small compressibility expansion was employed which requires the density 
fluctuations to be small.) There is an additional problem with the 
Smagorinsky model (16) even when it is restricted to flows with small density 
changes as pointed out by Bardina, et al.^ This model neglects the kinematic 
Leonard stress L^ and Cross stress given by 


J k£ 


v k v £ " 


v kV 
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( 20 ) 


(the total kinematic subgrid scale stress + R^ where 

s v k" v l' 1® kinematic Reynolds stress). For flows at low Mach 

numbers Mq < 0.1, density changes are extremely small and the difference 
between Favre filtered fields and ordinary filtered fields is negligible. 
Hence, at low Mach numbers, the model derived by Yoshizawa* would correspond 
to equation (12) with the scale similarity part (i.e., the first term on the 
right-hand side of (12)) neglected. In particular, it would correspond to the 
model 
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( 21 ) 


for the deviatoric part of the subgrid scale stress tensor. The 
computations conducted by Erlebacher, et al.^ indicate that such a model cor- 
relates poorly with the results of direct numerical simulations of compressi- 
ble isotropic turbulence at low Mach numbers. To be specific, correlations at 
the tensor level (i.e., correlations of exact versus modeled subgrid scale 
stress tensors) for Mq = 0.1 were only of the order of 30% as compared to 
the approximate value of 85% for the more complete model given by (12) (see 
Table 1 and 2). This improvement in the case of compressible flows is 
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comparable to the improvement obtained by the linear combination model for in- 
| compressible flows (see Bardina, et al.^). 

| A distinguishing feature of compressible turbulence models is the exist- 

ence of an isotropic subgrid scale stress which cannot simply be absorbed into 
the pressure term as is customarily done for incompressible flows. The iso- 
tropic part of Yoshizawa's model, in the limit of low Mach numbers for which 
the difference between Favre and ordinary filters are negligible, is given by 

i’m ‘ - T V i2lI ~ 6 M < 22 > 

s 

which differs from the Erlebacher, et al.^ model by the absence of the iso- 
tropic part of the Leonard and Cross stress contributions. However, this 
[ model correlates even more poorly than the deviatoric part ^ of the sub- 

grid scale stress tensor when compared with the results of direct numerical 
simulations of isotropic compressible turbulence. A scatterplot of the 
modeled (Eq. (22)) versus exact isotropic subgrid-scale stress j T kk^ Is 
shown in Figure 1 which was obtained from a 96^ direct simulation of iso- 
tropic compressible turbulence at Mq - 0.1. The correlation coefficient for 
| this case is 15%. This correlation improves to the 85% level when the Leonard 

! and Cross stresses are included.^ It must be kept in mind however, that di- 

\ rect numerical simulations (Erlebacher, et al.^) indicate that the isotropic 

Reynolds stress affects the correlations by at most a few percent for Mq < 

I 0.4. 

I 

In conclusion, it has been pointed out that the Smagorinsky-type subgrid 
scale model developed by Yoshizawa^ from the DIA for compressible large-eddy 
simulations is not applicable to turbulent flows with large density fluctua— 

| 

I 
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tions and correlates somewhat poorly with the results of direct numerical 
simulations of compressible isotropic turbulence at small Mach numbers* In 
order to overcome these limitations, we believe that models should be based on 
Favre filtered fields (where no restrictions are placed on the density fluctu- 
ations) and the Leonard and Cross stresses must be accounted for. The com- 
pressible subgrid scale models recently developed by the authors^ appear to be 
more general in this regard and are currently under study in a research pro- 
gram involving the large-eddy simulation of compressible turbulent shear 
flows • 
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M 0 = 

0.6 


L+C+R 
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Table 1. Comparison of correlation coefficients of the exact subgrid scale 
stress versus its model^ (given by Bq . (12)) ob- 

tained from a 96^ direct simulation of compressible isotropic turbulence 
(D = diagonal tensor level, OD = off-diagonal tensor level, V = vector level, 

S = scalar level, i -T<v t - C k* * ■* <v k''* + Vp- *kl ■ • 

Mq = Average Mach number over the computational domain) . 


M 0 


0.1 

0.4 

0.6 


D 

31 

31 

31 

ijR 

OD 

26 

26 

25 


V 

22 

22 

22 


s 

45 

45 

45 


Table 2. Comparison of correlation coefficients between the exact and the 
modeled deviatoric subgrid scale Reynolds stress (given by Eq. (21)) obtained 
from a 96 direct simulation of compressible isotropic turbulence 
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Figure 1. Scatterplot of the isotropic modeled 1 versus the exact subgrid 
scale stress obtained from a 96^ direct simulation of compressible 
isotropic turbulence at Mg * 0.1. 
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